† Corresponding author. E-mail:
Project supported by the National Innovation Fund of Chinese Academy of Sciences (Grant No. CXJJ-16M208), the Preeminent Youth Fund of Sichuan Province, China (Grant No. 2012JQ0012), and the Outstanding Youth Science Fund of Chinese Academy of Sciences.
We propose a new algorithm for wavefront sensing based on binary intensity modulation. The algorithm is based on the fact that a wavefront can be expended with a series of orthogonal and binary functions, the Walsh series. We use a spatial light modulator (SLM) to produce different binary-intensity-modulation patterns which are the simple linear transformation of the Walsh series. The optical fields under different binary-intensity-modulation patterns are detected with a photodiode. The relationships between the incident wavefront modulated with the patterns and their optical fields are built to determinate the coefficients of the Walsh series. More detailed and strict relationship equations are established with the algorithm by adding new modulation patterns according to the properties of the Walsh functions. An exact value can be acquired by solving the equations. Finally, with the help of phase unwrapping and smoothing, the wavefront can be reconstructed. The advantage of the algorithm is providing an analytical solution for the coefficients of the Walsh series to reconstruct the wavefront. The simulation experiments are presented and the effectiveness of the algorithm is demonstrated.
Wavefront sensing has been an indispensible part for adaptive optics (AO) systems. A typical adaptive optics system must have a way to sense the wavefront with enough spatial resolution and time resolution to apply a real-time correction.[1] The requirements of wavefront sensing in an adaptive optics system have a great difference from those of phase or figure determination in optical testing.[2] The wavefront sensing frequency can be on the order of hundreds of hertz. Meanwhile, with the wide applications of adaptive optics in a wide range from astronomical telescopes to microscopy, new challenges arise to meet the stringent requirements in applications, such as ground-based telescopes observing weak objects through the atmosphere.[3,4] Higher requirements on speed and sensitivity for wavefront sensing sensors are urgently necessary. However, most of the conventional wavefront sensing techniques use array detectors to acquire signals, such as Shack–Hartmann wavefront sensing,[5] curvature sensing,[6] and phase retrieval.[7] Two-dimensional detectors have significant limitations, especially in speed and sensitivity. They are also very expensive and hard to manufacture. In many cases, the array detectors working in some special conditions, such as weak light or high noise, often have low signal-to-noise ratios and cannot reconstruct the wavefront with enough accuracy. So another wavefront sensing approach measuring the wavefront with a single detector comes into consideration.[8,9] Compared with array detectors, single detectors ordinarily have higher performance and lower price. Due to the structure character, designs with single detectors can focus the light of the entire aperture and realize measurement without the light beam splitting and dividing. For this reason, the design of wavefront sensing with a single detector is likely to acquire a high signal-to-noise ratio under the condition of weak light. Besides, the signals detected with a single detector are one-dimensional, which means that the quantity of data is rather small. It provides a relative lower requirement on transmission and computing. Obviously, it is a potential heuristic direction sensing wavefront with a single detector.
There are two known papers realizing wavefront sensing using a single detector through the measurements of binary aberration modes.[10,11] The most outstanding contribution of them is providing a new way of wavefront sensing with a single detector. It can realize wavefront sensing with enough spatial resolution and speed which opens a great potentiality in adaptive optics system. Both of them expend the incident wavefront with a series of orthogonal and binary functions, the mode-field-Walsh (MOW) functions. The difference between them lies in that Feiling Wang manipulates the wavefront with binary phase patterns and Shuai Wang uses binary intensity patterns. A relationship equation between the coefficients of the MOW functions and their optical fields has been built. However, based on the iterative algorithm[3] or optimization strategy,[11] they can only obtain approximations of the coefficients by solving the equations. The low speed of the algorithms and the uncertain accuracy of the approximations limit the speed of the single detector as well as the method. So an algorithm that can provide exact coefficients with high speed is necessary and becomes the key point for the wavefront sensing method working well in real-time adaptive optics systems.
In this paper, we present a new algorithm which provides an analytic solution for the coefficients of the MOW functions. New selected modulation patterns derived from the linear transformation of the MOW functions are replenished for the algorithm. According to the modulation patterns, new relationships between the supplementary modulation patterns and their intensity signals (the light intensity within the pin hole) are built to determine an analytical solution for the coefficients of the MOW functions. A discrete wrapping wavefront whose phase is in the range of [
In order to realize wavefront sensing through intensity modulation, it is necessary to establish the connection between intensity and phase. It should be noted that no light is a special situation in intensity modulation which means that the amplitude and phase of light disappear at the same time. Meanwhile, when the amplitude of the modulated light is not zero, the phase of the modulated light keeps the same as that of the incident light. So the intensity modulation produces two kinds of special statements, “light on” meaning that the phase is chosen and “light off” meaning the phase is ignored.
The set of Walsh functions is a series of complete and orthogonal binary functions that take on two values, −1 and +1. The MOW functions are derived from the two-dimensional Walsh functions in the polar coordinates. The Walsh functions are formed by dividing the azimuth angle and the radius. The division is used to produce an orthogonal and complete series in the circular aperture. figure
The Walsh functions have an important property
The optical arrangement of the wavefront sensing is shown in fig.
According to the Fourier diffraction theory, the signal measured at the ideal focal point by the photodiode is given by[8]
We expand the aberration function of the incident wavefront with the MOW functions in the following form:
In order to acquire the coefficients of the MOW functions, we need to build connections between the modulation patterns and the MOW functions. The intensity SLM modulator can only influence the amplitude of the incident light directly. As discussed above, the intensity modulation produces two kinds of special statements, light on and light off. So a set of relationships between the binary-intensity-modulation patterns and the MOW functions can be built by a simple linear transformation[11]
The wavefront modulated with the binary-intensity-modulation patterns can be described as
As discussed above, the phase of the reconstructed wavefront is
Because
We choose the Hadamard order to produce the MOW functions in the wavefront sensing method. According to Eq. (
In Eq. (
By now, the coefficients
It should be noticed that the value of the wrapped wavefront phase is limited in the range of
In order to recover the wavefront, the sampling rate should satisfy the Nyquist sampling theory. So phase unwrapping techniques should start from the fact that the phase differences of the neighboring pixels are less than π.[12] And the unwrapping operation is relatively easy to implement when the sampling rate is higher. The unwrapped relative phase can be written as
We produce the incident wavefront with the first 35 Zernike polynomials depending on the power spectrum of the Kolmogorov model for atmospheric turbulence.[14] figure
The wavefront sensing frequency is important in an AO system. The wavefront sensing algorithm needs the intensity signals of the wavefront modulated with a series of modulation patterns. The number of the patterns means the time cost. So in order to increase the wavefront sensing frequency, the number of the MOW functions should be limited. We reconstruct the incident wavefront with the first 64 Walsh functions. The reconstruction result is shown in fig.
Compared with the wavefront reconstructed with the first 256 Walsh functions, the wavefront reconstructed with the first 64 Walsh functions has a large residual. So the number of the Walsh functions used has a great impact on the reconstruction precision. Besides, some special Zernike orders cannot be reconstructed with the first 64 Walsh functions, such as the 13th and 24th orders. The wavefront generated with the first 35 Zernike polynomials cannot be reconstructed accurately with the first 64 Walsh functions.
The wavefront sensing algorithm can work well over different shapes. In fact, it can work well with different shapes of the pupil under the condition that the pupil can be divided into 2n (the index n can be any integer) units and the square of each unit is equal. We make a simulation to illustrate it. Figure
Meanwhile, we reconstruct the wavefront with the stochastic parallel gradient descent (SPGD) algorithm. Firstly, we can detect the light intensities, I, of the light beam under different modulation patterns. Then in the parallel stochastic optimization, we assume that the initial values of the coefficients of the first 256 MOW functions are zeros. Next a two-side random ensemble small perturbation
As shown in fig.
Benefitting from the high sensing speed and sensitivity of a single detector, sufficient signals can be acquired for wavefront sensing in an adaptive optics system. It provides a hardware foundation for the wavefront sensing method based on binary aberration modes applied in a real-time AO system. However, the coefficients of the MOW functions cannot be acquired with enough speed and accuracy by the iterative algorithm or optimization strategy. It restricts the use of the wavefront sensing method in adaptive optics.
In this paper, a new algorithm which can provide an analytical solution for the wavefront sensing method based on binary aberration modes is proposed. New modulation patterns are produced to build detailed and strict relationships. The new modulation patterns can help to acquire the coefficients of the Walsh functions. A total of
As the algorithm can realize the reconstruction of the wavefront phase accurately, the wavefront sensing method could also be used in optical testing. A higher accuracy and a larger dynamic range can be achieved with more modulation modes. Besides, the method can sense a phase stage. So it may have a potential application in the reconstruction of the wavefront containing step jumps.
The algorithm acquires an analytical solution by adding new modulation patterns, which would require a relatively long time to achieve one measurement. This produces a limitation for its applications in a real-time AO system. It makes the wavefront sensing frequency of the method lower compared with that of the traditional method, such as the Shack–Hartmann wavefront sensing method. However, the frequency can be improved by compressing the number of modulation patterns. According to the Fourier diffraction theory, the sign of wavefront phase φ cannot be determined with the wavefront sensing method. Corresponding methods on how to determine the sign of the wavefront are necessary for the wavefront sensing method showing its potentiality in adaptive optics systems. The wavefront sensing algorithm still has some open questions. Further study is required to complete it.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] |